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^ ABSTRACT 

Q Distant galaxy clusters provide important tests of the growth of large scale structure 

^ in addition to highlighting the process of galaxy evolution in a consistently defined 

C/^ environment at large look back time. We present a sample of 22 distant [z > 0.8) 

^ galaxy clusters and cluster candidates selected from the 9 deg^ footprint of the over- 

' ' lapping X-ray Multi Mirror (XMM) Large Scale Structure (LSS), CFHTLS Wide 

and Spitzer SWIRE surveys. Clusters are selected as extended X-ray sources with an 

^ accompanying overdensity of galaxies displaying optical to mid-infrared photometry 

1/^ consistent with z > 0.8. Nine clusters have confirmed spectroscopic redshifts in the 

00 interval 0.8 < z < 1.2, four of which are presented here for the first time. A further 

11 candidate clusters have between 8 and 10 band photometric redshifts in the in- 
terval 0.8 < z < 2.2, while the remaining two candidates do not have information in 
sufficient wavebands to generate a reliable photometric redshift. All of the candidate 

i-H clusters reported in this paper are presented for the first time. Those confirmed and 

01 candidate clusters with available near infrared photometry display evidence for a red 

sequence galaxy population, determined either individually or via a stacking analysis, 

^ whose colour is consistent with the expectation of an old, coeval stellar population 

• i-H observed at the cluster redshift. We further note that the sample displays a large 

range of red fraction values indicating that the clusters may be at different stages 
^ of red sequence assembly. We compare the observed X-ray emission to the ffux ex- 

pected from a suite of model clusters and find that the sample displays an effective 
mass limit M200 ~ 1 x lO^^M© with all clusters displaying masses consistent with 
^200 < 5 X 1O^^M0. This XMM distant cluster study represents a complete sample of 
X-ray selected z > 0.8 clusters. We discuss the importance of this sample to investi- 
gate the abundance of high redshift clusters and to provide a relatively unbiased view 
of distant cluster galaxy populations. 
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Observations of galaxy clusters provide crucial insight into 
the development of structure in the Universe, from the 
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growth of clusters themselves, to the evolution of their mem- 
ber galaxies. Furthermore, cluster studies yield important 
constraints on cosmological models through tests of the 



direct evidence is emerging of both red sequence truncation 
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growth of structure (e.g. Vikhlinin et al. 
aL]|2011 ). The greatest constraining power on cosmological 



parameters and on the co-evolution of galaxies and clus- 
ters requires large, well- controlled samples of clusters out 
to ^ > 1. To date, a number of techniques have been suc- 
cessfully applied in order to generate such samples of clus- 
ters at 2; < 1. These include approaches based upon de- 
tecting galaxy overdensities (in a combination of apparent 
colour and sky position; e.g. [Postman et al.|p"996| |Gladders| 
& Yee 2005), extended X -ray sources (e.g.|Gioiaet al.|1990[ 



Bohr inger et al. 2000! [Burenin et al.|| 2007| [Mehrtens et aT] 
2012| ) and the Sunyaev-Zel'dovich decrement observed to- 



ward the Cosmic Microwave Background (e.g. Menanteau 
et al.||2010||Reichardt et al.||2012| ) 



Although the systematic estimation of the cluster num- 
ber density above z > 1 is an important issue the search 
for high redshift clusters is made difficult by the faintness 
of the cluster signal, e.g. small galaxy over-density in opti- 
cal and near-infrared (NIR) imaging, weak X-ray emission 
whose extent is difficult to assess, etc. However, should X-ray 
imaging observations of sufficient depth and spatial resolu- 
tion be executed, the detection of high redshift clusters via 
spatially-extended emission is advantageous as it provides 
clear evidence of hot gas confined within a gravitational po- 
tential well. Care must be taken though to assess the extent 
to which faint X-ray active galactic nuclei (AGN) may mimic 
or modify the significance and spatial extent of cluster emis- 
sion and its spectral form (e.g. Branchesi et al.|2007 Pierre 
et al.||2012t 



Systematic searches for high redshift galaxy clusters in 
X-rays are currently being conducted via dedicated XMM 
surveys (XDCP|Fassbender et al.|2011| XMM-LS S [Pierre et 
|al.|2007|[Fremer et al.|2006| [Andreon et al.|20"05l , this work; 
XCS [Romer et al.|2001| [Stanford et al.|2006| ). These surveys 
employ quantitative algorithms to identify extended sources 
and are characterised by accurate selection functions and 
a clearly defined relationship between cluster observables 
(such as X-ray luminosity) and total cluster mass ("Reichert 
jet al.|20lT ). Galaxy clusters at 2; > 1 may also be detected by 
the extension of successful "red sequence" searches to near- 
and mid-infrared (MIR) wavebands in order to detected the 
redshifted emission from distant galaxies (e.g. [Muzzin et al 



2009 ) . An associated technique employed by the IRAC Shal- 



low Cluster Survey (Eisenhardt et al. 2008) identifes clusters 
to z < 1.5 as stellar mass overdensities in multi-band pho- 
tometric redshift slices. Each MIR imaging technique has 
proven very successful at identifying large numbers of can- 
didate and confirmed clusters. 

Galaxy clusters ai z < 1 have been employed exten- 
sively to study their member galaxy populations and indi- 
cate that they are composed of uniformly old stellar pop- 
ulations where the bulk of their stars formed at 2; = 3 or 
greater ( [Jaffe et al.[|2011| . In addition, low redshift clusters 
display strong population trends such as the morphology- 
density relation (e.g. Dressier et al.|19"97 ). Such relations re- 
flect the dominance of bright, red, bulge-dominated galaxies 
in cluster cores. Observing clusters at high redshift provides 
an opportunity to approach the epoch when the progenitors 
of low-redshift galaxies where assembled. Indeed impressive 



(Rudnick et al. 2012) and merger driven galaxy assembly 



(Lotz et al. 2012) in a MIR selected (yet X-ray detected) 
cluster at z — 1.6. The evolution of the brightest cluster 
galaxies (BCGs) in distant clusters presents a more com- 
plex picture: Lidman et al.| (^12) report that the stellar 
mass of BCGs in 160 clusters spanning 0.03 < z < 1.63 
grow steadily with decreasing redshift in a manner consis- 
tent with a semi-analytic model. In potential contrast to 
such evolution in high redshift cluster galaxies, [Stott et al.j 
(2010 ) report that the stellar mass contained in the BCGs 
of a sample of 20, 0.8 < z < 1.5 clusters has changed lit- 
tle between the epoch of observation and the present day. 
Such trends observed in heterogeneously assembled samples 
will be better understood by performing similarly extensive 
analyses upon a complete sample of high redshift clusters 
selected employing a single method (e.g. Fassbender et al 



2011). 



The XMM-LSS survey is well placed to contribute to 
this investigation: covering 11 deg^ with X-ray imaging to a 
depth of ^ 1 X 10~^^ erg s~^ cm~^ for extended sources in 
the [0.5-2]keV waveband and accompanied by optical and 
MIR photometry. The XMM-LSS project has previously 
demonstrated the ability to identify clusters to z = 1.2 (Br 
mer et al. 2006) and has published z < 1 cluster number 



counts selected according to a clear, quantitative selection 
function ( [Pacaud et al.[[2007 ). 



This paper presents the XMM-LSS distant cluster sam- 
ple. These are defined to be extended X-ray sources at 
2; > 0.8 and consist of spectroscopically confirmed clusters 
together with a number of candidate clusters supported by 
a detailed photometric redshift analysis. The distant cluster 
sample has been identified from the full sample of extended 
X-ray sources within a 9 deg^ sub-area of the XMM-LSS 
survey and in this sense it represents a complete sample 
of X-ray selected distant clusters. In particular, the meth- 
ods employed to determine whether a given extended X-ray 
source is or is not a distant cluster are selected to minimise 
any potential bias such as the presence of a strong red se- 
quence. In this sense the sample should be as complete as 
possible and should provide an unbiased perspective of the 
galaxy populations in distant X-ray clusters. 

The structure of the paper is as follows. Section 2 de- 
scribes the construction of the distant cluster sample (con- 
taining both spectroscopically confirmed and candidate sys- 
tems) and the multi- wavelength data sets used to define it in 
addition to presenting images of the sources in the sample. 
Section 3 describes the approaches taken to explore which of 
the candidate systems have clear photometric evidence for 
being genuine distant clusters. Section 4 discusses the results 
of applying these approaches and explores the diversity of 
properties shown by the likely high redshift systems. 

This paper employs a Friedmann-Robertson- Walker 
cosmological model described by the parameters Qm = 0.3, 
= 0.7, Ho = 70kms~^Mpc~^. An angular scale of 1' cor- 
responds to a transverse physical scale of 480 kpc, 508 kpc 
and 502 kpc at redshifts z = 1, 1.5 and 2 respectively. All 
photometry is quoted in the AB system. 
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2 THE XMM-LSS SURVEY SAMPLE 

The XMM-LSS survey currently covers approximately 11 
deg^ and is described in Chiappetti et al.| (|2012| and Clerc 
et al. (2013). Galaxy clusters are detected as extended X- 
ray sources and are classified as either CI or C2 on the 



bas is of their surface brightness characteristics (Pacaud et 

' T 



-14 

erg s cm 



aL[2006). The effective flux limit is - 1 x 10" 
for extended sources. 

Visual inspection of the X-ray images of individual sys- 
tems along with their optical and NIR images confirm that 
the CI class represents an uncontaminated sample of ex- 
tended X-ray sources (mainly clusters but with some detec- 
tions of X-ray halos in very low redshift galaxies, |Pacaud| 
et al. 12006). The C2 class displays a contamination rate of 



30-50%, with the main sources of contamination being mis- 
classified point sources and artefacts on the X-ray image. 
These contaminants are typically removed by visual inspec- 
tion of the X-ray image prior to further analysis. 

The 11 deg^ XMM-LSS sample contains 50 CI and 60 
C2 sources, of which 44 CI and 27 C2 sources have confirmed 
redshifts from optical spectroscopy. The redshift distribution 
of confirmed sources ranges over 0.05 < z < 1.22. The lower 
spectroscopic confirmation rate for the C2 sources arises due 
to a) the lower priority placed on follow-up of such sources 
compared to CI sources and b) the increased difficulty of 
following up fainter, lower quality detections. 

The analysis used to generate the XMM-LSS distant 
cluster sample is based upon a 9 deg^ subregion of the 
XMM-LSS field. This region represents the common foot- 
prints of the XMM-LSS, Canada France Hawaii Telescope 
Legacy Survey (CFHTLS) Wl and Spitzer SWIRE 3.6 and 
4.5^Lm surveys and contains 88 C1+C2 sources (of which 
55 have spectroscopic redshifts). NIR imaging drawn from a 
variety of sources (i.e. the UKIDSS and WIRDS surveys in 
addition to individual CFHT/WIRCam and VLT/HAWKI 
images) exists for many of the spectroscopically confirmed 
and candidate distant systems, with the available bands and 
depth varying on a source- by source basis. The principal 
data sets and their processing are described below. 



2.1 X-ray photometry 

The characterisation and measurement of extended sources 
in the XMM-LSS survey is described in detail in |Pacaud| 
|et al.| (^2006 ) and [Adami et al.| ([2011| and we summarise 
the important features of the analysis here. Sources are de- 
tected above a specified pixel threshold by appliying the 
SExtractor routine to a multi-scale wavelet reconstruction 
of the XMM science image. Individual sources are char- 
acterised as either extended or point-like on the basis of 
the likelihood values of appropriate models applied to each 
source. Taking the example of an extended source, nearby 
point-sources identified individually by the detection algo- 
rithm are masked when performing the source extended 
fit using the SExtractor segmentation map ( [Pacaud et al 



[200 6 ). Therefore the extent likelihood used for C1/C2 clas- 
sification is almost free from contamination by bright point 
sources. Faint point-sources contaminating the extended 
source emission and not deblended by the algorithm may 
be present, particularly in regions close to the cluster cen- 
tre. Accounting for them is challenging given the faint ness 



of our objects and their compactness relative to XMM point 
spread function. 

We compute X-ray fluxes for the extended sources asso- 
ciated with spectroscopically confirmed and candidate clus- 
ters in the [0.5-2] keV band employing the procedure out- 
lined in Adami et al. (2011 ). The method applies a curve-of- 



growth analysis to the X-ray count rate data which confers 
the advantage of being free of any profile assumptions ap- 
plied to the extended X-ray source. We estimate that the 
application of a finite size aperture used for these measure- 
ments recovers 80-90% of the total count rate of the cluster 
and we note that this bias is lower than the statistical error 
of our measurements. The flux measurement step allows a 
further check for additional blended emission. Point sources 
lying close to but off the cluster central region are either 
identified by the detection algorithm or flagged visually. In 
both cases their contribution to the total flux is excluded 
from the extraction region, and the missing area is accounted 
for by assuming a circularly symmetric flux profile. The 
fluxes were obtained assuming a fixed conversion factor of 
9 X 10^^(ergs~^cm~^)/(cnts s~^). This value was calculated 
using Xspec from an APEC emission model with the follow- 
ing parameters: z = 1, T = 4 keV, Nh = 2.6 x 10^°cm^, 
Ab = 0.3 (note that the conversion factor changes by less 
than 5% for z = 1.5 and T = 4 keV). Bolometric luminosities 
were calculated with Xspec employing the measured fluxes, 
the cluster redshift (quoted in Table n!| and the fitted cl us- 
ter temperature (B remer et al.|20 06^; Paca ud et al.|2007| or 
assuming T — 4 keV if no temperature was available. Flux 
and luminosity values for each cluster are listed in Table [l] 

2.2 Multi-wavelength photometric data 

As noted above, most of the data used in this analysis are 
taken from existing large area surveys which cover part or 
ah of the XMM-LSS region. Descriptions of the CFHTLS 
and Spitzer SWIRE data used in the paper are described 
in Gwyn (2012) and Chiappetti et al. (2012) respectively. 
Large area NIR survey data provided J, H and/or iT-band 
imaging for a number of the CI and C2 z > 0.8 candi- 
dates. These data came from either the UKIDSS Deep Ex- 
tragalactic Survey (DXS DR8; see Lawrence et al.||2007 ) or 
the WIRCam Deep Survey (WIRDS; [Bielby et al.|[2010 ). 
Typical depths (5cr) for these data are J = 23.4, K = 22.9 
(DXS) and J = 24.7, if = 24.7, i^^s = 24.7 (WIRDS). NIR 
imaging data for the remaining sources in the distant cluster 
sample were obtained with principal investigator programs 
using VLT/HAWKI and CFHT/WIRCam and are described 
below. Details of the available multi- wavelength photometry 
for each z > 0.8 system, both candidate and confirmed, are 
given in Table [l] 
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Table 1. Distant clusters in the XMM-LSS survey. These sources represent all C1+C2 sources (Section 2) within the approximately 9 deg^ XMM-LSS/CHFTLS/SWIRE footprint that 
cannot be unambiguously associated with a 2; < 0.8 cluster or galaxy. Spectroscopically confirmed distant clusters are displayed in the upper section of the table and have their redshift 
values indicated. Each is also described by an official XLSSC identifier. Candidate distant clusters are displayed in the middle section of the table and their photometric redshift values 
(where computed) are indicated with accompanying uncertainty. Unknown sources are displayed in the lower section of the table and do not have a redshift value. Each cluster is further 
identified using a sequential number (column 1) to simplify their discussion in this paper. In subseqent publications that refer to these sources, we recommend using either the cluster 
names or the XLSSC numbers which are the only lAU validated labelings. The procedures used to compute the X-ray fiux and luminosity measures are described in Section 2.1. 
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Eight of the CI and C2 X-ray sources that were can- 
didate or confirmed z > 0.8 clusters were observed us- 
ing HAWK-I on the VLT through ESQ programme 084.A- 
0740(A). The HAWK-I camera consists of four HAWAII 
2RG 2048 X 2048 pixel detectors. The four detectors image 
an area of 7.^5 x 7^5 with a pixel scale of 0Vl06 pixeP^ (the 
cross shaped gap between the four detectors is 15" wide). 
Each candidate was observed using the YJKs filters with 
exposure times of 1800 seconds, 2530 seconds, and 3600 sec- 
onds respectively and a suitable offset was applied to each 
set of observations to place the measured X-ray centroid 
of each cluster candidate in the centre of a single detector. 
Data were obtained in service mode during October 2009 to 
January 2010. 

Reduction was carried our using standard procedures 
from version 1.4.2 of the ESO pipeline within the esorex 
environment. The pipeline routines corrected the data for 
the presence of bad pixels, dark current, fiat field variations, 
two-stage sky subtraction with object masking, distortion 
correction and co- addition with pixel rejection. 

Photometric calibration was a two-step process. Firstly, 
standard star observations were used to place all four detec- 
tors on a common photometric scale. Then the J and Ks- 
band images were placed on the 2MASS scale by matching 
stars with J < 15.8 and Kg < 14.3 to their counterparts 
in the HAWK-I images. Between 3 and 20 such stars were 
present in each image and this procedure resulted in zero- 
points accurate to better than 0.1 mags in all cases. The of- 
ficial HAWK-I y-band zeropoint is yet to be made available. 



In its absence, the procedures used in Hickey et al. (2010) 



were used. Using the HAWK-I J-band photometry and the 
CFHTLS 2;-band, a pseudo F-band magnitude was con- 
structed from their fiux average for each source. Using only 
those sources with z — J ^ 0, the HAWK-I Y zero point was 
adjusted to match the pseudo-F photometry. The typical 
sensitivity of these data are Y = 24.3, J = 24.0, Kg = 23.5 
(5cr) within a 3^^ diameter aperture. 

A further 4 candidate C2 distant clusters were observed 
with CFHT/WIRCam in December- January 2011-12. The 
WIRCam camera consists of four HAWAII2-RG detectors, 
each containing 2048 x 2048 pixels. The four detectors image 
an area of 20^ x 20^ with a pixel scale of 0^3 pixel" ^ (the cross 
shaped gap between the four detectors is 45^^ wide). Each 
candidate was observed using the J Kg filters with exposure 
times of 5428 to 8614 seconds, and 3675 seconds respectively. 
A suitable offset was applied to each set of observations to 
place the measured X-ray centroid of each cluster candidate 
in the centre of a single detector. Photometric calibration 
followed the method outlined above and generated typical 
sensitivity values J = 23.0,^^5 = 22.5 (5cr) within a 3^^ 
diameter aperture. 



each cas^ Photometry was computed within an aperture 
based upon the Kron (1980) radius. The image quality of 
the HAWK-I and WIRCam data are well matched to that 
of the CFHTLS optical data, with typical stellar FWHM 
values of 0'.'6 - 0'.'7. 

Finally, we match the optical-NIR catalogues for each 
field to sources detected at 3.6|J,m and 4.5|J,m in the SWIRE 
catalogue. Sources are matched with a 1^.^5 tolerance and 
Spitzer photometry is quoted within a 2" radius circular 
aperture with an additive offset applied to correct to a 
pseudo-total aperture. We note that this approach may in- 
troduce a small additive offset between the matched optical- 
NIR and the Sptizer photometry and we attempt to deter- 
mine and correct for any such zero-point offsets between 
wavebands within the photometric redshift analysis. 



2.4 Spectroscopic observations 

Nine galaxy clusters presented in this distant cluster sample 
have confirmed spectroscopic redshifts. Five have been pub- 
lished previously and appropriate references are provided in 
Table [l] Spectroscopic data for four of the clusters are pre- 
sented here for the first time and we decribe the observation, 
reduction and analysis of the data below. 

Details of the spectroscopic observations obtained for 
clusters 01, 02, 04 and 06 are presented in Table [2] Slit tar- 
gets for each cluster consisted of galaxies located within the 
X-ray emitting isophotes with photometric redshifts consis- 
tent with being cluster members. Further slits were placed 
upon moderately bright stars and galaxies in order to pro- 
vide identifiable reference spectra. The data were processed 
using standard techniques in the IRAlJ^ environment which 
included procedures to extract one dimensional spectra, ap- 
ply a wavelength solution based upon reference HeAr spectra 
(or sky features in the case of GMOS-S spectra) and to cor- 
rect for varying spectrograph response using observations of 
a spectrophotometric standard star. The spectral resolution 
of each data set was determined by measuring the full- width 
at half- maximum (FWHM) of unresolved arc lines. 

Reduced spectra were inspected visually to provide an 
initial estimate of galaxy redshifts based upon the identifi- 
cation of prominent features. Individual spectra were then 
cross-correlated with a representative early-type galaxy 
template (e.g. Kinney et al.|1996 ) employing routines based 
upon [Tbnry Davis (1979). Individual cluster members 
were selected by identifying visually an initial cluster red- 



shift Zpeak in the redshift histogram for each field. Cluster 
galaxies were then selected to display Zpeak ±0.03 and indi- 
vidual redshift values are displayed in Table |3] The cluster 
redshifts presented in Table ^ were then computed as the 



2.3 Catalogue construction 

Optical and NIR imaging data for each X-ray source were 
placed on a common pixel scale using the swarp v2 . 17 . 1 
software package (Bertin et al. 2002). Source extraction and 
photometry were then performed using SExtractor v2 . 5 . 
( Bertin Arnouts]|1996| ) employed in dual- image mode us- 
ing the Ks- or K-band image as the detection image in 



^ Note that for brevity in the following text we refer to Kg 
photometry for all relevant sources as this is the K-band fil- 
ter predominantly used and as the colour term between the 
UKIDSS K filter and the Kg filters used for CFHT/WIRCam 
and VLT/HAWK-I is small for the z > 0.8 galaxies of interest in 
this paper. 

^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Table 2. Details of the spectroscopic observations performed on four clusters in the sample. 



Lyluster 


Observing date 


Telescope / Spectrograph 


Grism+Filter 


Wavelength 


spectral 


Exposure 










coverage (A) 


resolution (A) 


time (s) 


01 


November 2003 


VLT/FORS2 


600RI+GG435 


5000-8500 


7 


3600 


02 


November 2003 


VLT/FORS2 


600RI+GG435 


5000-8500 


7 


3600 


04 


November/December 2006 


Gemini/GMOS-S 


R400+OG515 


6000-10000 


8 


17500 


06 


November 2010 


VLT/FORS2 


300I+OG590 


6000-10000 


20 


16200 



Table 3. Spectroscopic redshifts for individual galaxies in each 
cluster. 



ID 


R.A. (J2000) 


Dec. (J2000) 


redshift 


01 


35.9829 


-3.3481 


0.799 


01 


35.9870 


-3.3455 


0.801 


01 


36.0002 


-3.3564 


0.800 


01 


36.0026 


-3.4265 


0.799 


01 


36.0042 


-3.4272 


0.807 


01 


36.0003 


-3.4263 


0.803 


01 


36.0373 


-3.3946 


0.801 


01 


36.0108 


-3.4389 


0.803 


02 


35.6420 


-4.9655 


0.832 


02 


35.6430 


-4.9689 


0.825 


02 


35.6397 


-4.9720 


0.842 


02 


35.6401 


-4.9587 


0.832 


02 


35.6414 


-4.9775 


0.843 


02 


35.6345 


-4.9761 


0.833 


02 


35.6832 


-5.0105 


0.831 


02 


35.6714 


-5.0079 


0.822 


02 


35.6121 


-4.9991 


0.861 


02 


35.6657 


-5.0137 


0.844 


02 


35.6114 


-4.9987 


0.862 


02 


35.6413 


-4.9652 


0.834 


02 


35.6824 


-5.0102 


0.832 


04 


33.8869 


-3.7395 


1.000 


04 


33.8826 


-3.7344 


1.008 


04 


33.8801 


-3.7197 


1.004 


04 


33.8771 


-3.7355 


0.999 


04 


33.8593 


-3.7478 


1.010 


04 


33.8571 


-3.7290 


1.002 


04 


33.8535 


-3.7269 


1.006 


04 


33.8494 


-3.7268 


1.003 


04 


33.8461 


-3.7275 


0.989 


04 


33.8435 


-3.7203 


0.999 


04 


33.8412 


-3.7321 


1.002 


04 


33.8352 


-3.7298 


1.011 


04 


33.8231 


-3.7332 


1.000 


04 


33.8070 


-3.7649 


0.999 


06 


33.7371 


-3.5019 


1.033 


06 


33.7406 


-3.4999 


1.032 


06 


33.7439 


-3.5019 


1.029 


06 


33.7460 


-3.5062 


1.033 



mean redshift of all members located within l' of the clus- 
ter X-ray centroid. 



3 IDENTIFYING DISTANT CLUSTER 
CANDIDATES 

The aim of this paper is to present a complete sample of 
extended X-ray sources with redshifts z > 0.8. The moti- 
vation for generating a complete sample essentially centres 
upon the assertion that, only by generating a complete sam- 
ple of X-ray selected ^ > 0.8 clusters can one a) compare 
the abundance of distant clusters to a cosmological model 
prediction in quantitative manner and b) discuss the range 
of galaxy properties exhibited by the member population in 
a relatively unbiased way. Constructing a complete sample 
ultimately requires providing a robust explanation of the 
nature of every extended X-ray source in the sample area. 
Though it is the aim of the XMM-LSS survey to confirm 
spectroscopically all extended CI and C2 sources, this re- 
mains an observationally challenging prospect at this stage. 
The methodology of this paper is therefore to employ the 
available photometric data to direct the detailed follow-up 
of individual clusters (e.g. deep NIR imaging and spectro- 
scopic observations) towards extended sources that show 
compelling evidence for being a bona-fide cluster above some 
redshift threshold. 

The threshold of z > 0.8 applied in this paper was 
defined in response to largely practical considerations, e.g. 
given the depth of available optical data, experience indi- 
cated that all clusters at z < 0.8 could be recognised with 
little ambiguity. A threshold of z > 0.8 (as opposed to a 
larger value) generated a final sample containing approxi- 
mately 20 distant clusters which was deemed large enough 
that evolution of cluster galaxy properties such as red se- 
quence colour and population mix could be traced within 
the sample rather than via reference to other cluster sam- 
ples selected possibly using alternative techniques. Finally a 
threshold of z > 0.8 corresponds to a look back time of ap- 
proximately 7 Gyr and accords with the redshift definition 
applied to distant clusters in the literature. 

We summarise below the steps employed to evaluate the 
robustness (or otherwise) of the evidence pointing to a high 
redshift identification for each candidate: 

• Classify visually all C1/C2 sources based on rz3.6\im 
and 3.6^im4.5^im images. Identify potentially distant sys- 
tems and reject misclassified point sources. 

• Compute the surface density of galaxies along the line 
of sight to each source that satisfy photometric selection cri- 
teria appropriate for identifying z > 0.8 galaxies. Compare 
to the visual classification results as a check. 

• Compute the colour of galaxy overdensities along the 
line of sight to each source and compare to the expected 
colour of model galaxies as a test of the distant cluster hy- 
pothesis. 
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• Compute photometric redshifts for galaxies along the 
line of sight to distant cluster candidates with available 
multi- wavelength photometry. Identify photometric redshift 
peaks spatially associated with the location of the extended 
X-ray source. 

• Obtain where possible spectroscopic observations of dis- 
tant cluster candidates with Zphot > 0.8 with the aim of con- 
firming at least 3 concordant redshifts within the sky area 
giving rise to the extended X-ray emission. 




3.1 Visual classification 

The first step is straightforward and is carried out on all CI 
and C2 sources in the 9 deg^ XMM-LSS/CFHTLS/SWIRE 
field. The CFHTLS and SWIRE fields containing each CI 
and C2 X-ray source were inspected visually, in order to 
identify any obvious clustering at or close to the X-ray po- 
sition. Images in individual bands and pseudo-true colour 
images (r / z/3.6\im and r/3.6^Lm/4.5|J,m) were used in this 
process. This step was performed by up to six people and 
the final classification was subject to the decision of two 
moderators. 

Visual inspection is relatively rapid to perform and pro- 
vides useful information on the broad nature of each X-ray 
source e.g. bright, clustered galaxies consistent with a low 
redshift cluster, a misclassified point source, or a significant, 
extended X-ray source with at best a grouping of faint, 
red galaxies consistent with being a high-redshift cluster. 
However, visual classification is likely to provide only an 
imprecise estimate of the redshift of each system. In prac- 
tice, though all 88 C1+C2 sources were inspected visually, 
the classification efforts were focussed on the 33/88 C1+C2 
sources lacking a secure spectroscopic identification. 

Of these 33 sources, 12 displayed convincing evidence 
for being a cluster at a redshift sufficiently below z = 0.8 
to be a secure classification (Figure [l] provides an image of 
a typical system). The remaining sources displayed either 
a) weak or absent optical emission with an identifiable clus- 
tering of 3.6^im sources, b) a clear case of a misclassified 
point source, or c) a weak, yet potentially extended X-ray 
source with no identifiable overdensity of galaxies at 3.6|am. 
Of these classes of objects, those in group (a) are retained 
as candidate z > 0.8 clusters (we note at this stage that 
the visual classification may include sources at z < 0.8 but 
the philosophy at this stage is that any cut avoid being con- 
servative). Sources in group (c) could either represent false 
detections resulting from the X-ray pipeline or potentially 
very distant X-ray clusters with little evidence of galaxy 
clustering at any wavelength. In either case, the prospects 
for spectroscopic confirmation of such systems are excep- 
tionally poor and - as they represent at best 1 or 2 systems 
out of a sample of 88 - are not retained as candidate distant 
clusters at this stage. 

This simple sifting can work to at least z — 1.2 as the 



highest redshift XMM-LSS cluster XLSSC 046 (Bremer et 



al. 2006) is straightforwardly selected in this way with an 



obvious compact overdensity of 3.6 ^im galaxies at the X-ray 
position with no bright optical counterparts (see Figure |2]). 
However, at high redshifts whether a cluster is straightfor- 
wardly identifiable depends on several factors affecting its 
surface density contrast against the background and fore- 
ground galaxies. If the foreground/background density in 




Figure 1. An example image of one of the 12 extended X-ray 
sources unambiguously associated with a 2; < 0.8 cluster on the 
basis of visual inspection. The image is composed of rz?>.Q\im data 
and is 2' on a side with standard astronomical orientation. The 
white contours indicate X-ray emission. 




mm 



Figure 2. A comparison of the appearance of the confirmed 
z = 1.2 cluster XLSSC 046 in r-band (left panel) and the 3.6|J.m 
waveband (right panel) . Each image is approximately 2' on a side 
with standard astronomical orientation. The white contours indi- 
cate X-ray emission. 



the z — 1.2 cluster field had been higher and the spatial dis- 
tribution of red cluster galaxies less compact (spread over a 
1' radius region rather than the 15^^ radius of the detected 
overdensity), the system would have been harder to discern 
in this way. 

Two exceptions to this process are clusters 06 and 19. 
Cluster 06 lies at the very edge of the SWIRE footprint and 
has data available at 3.6^Lm but not 4.5^Lm. Cluster candi- 
date 19 lies just outside the SWIRE footprint yet had previ- 
ously been flagged as a high redshift candidate on the basis 
of the extended X-ray image and faint, i-band detection of 
the candidate BCG. We include it in the following discus- 
sion and in Table ^ yet do not include it in the number of 
systems quoted in the 2; > 0.8 area limited sample. 

The visual classification supported by the existing spec- 
troscopic observations populates each cluster class as fol- 
lows: spectroscopically confirmed or candidate clusters at 
z < 0.8, 57 ("z < 0.8 clusters" hereafter); spectroscopi- 
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cally confirmed clusters at z > 0.8, 9 (see Figure [s] "con- 
firmed clusters" hereafter); candidate clusters at z > 0.8, 14 
(see Figure |4] "candidate clusters" hereafter); point-like or 
marginal sources with no clear identification, 8, (see Figure 
[sj "unknown extended sources" hereafter). 

3.2 The surface density of faint, red galaxies 

Following the visual inspection, the colours of galaxies in 
all fields, whether or not identified as clustered from initial 
inspection, are examined in order to identify any clustering 
in both position and colour. These colours, given data of 
suitable depth, can provide an indication of the redshift of 
the system. 

The photometric analysis employs the r — 3.6 ^im and 
3.6^Lm — 4.5\im colours derived from the CFHTLS and 
SPITZER/IRAC data. These colours have been proven to 
select for distant cluster galaxy populations (see Papovich 
2008||Muzzin et al.|200'8t . With data of sufficient depth, the 



use of the 3.6|J.m — 4.5|J.m colour is particularly powerful as 
it is a strong function of redshift between 0.5 < z < 1.5 for 
galaxies with a wide range of star formation histories and 
the scatter in colour between histories is relatively small 
(see Papovicli||2008| . For all clusters, the surface density of 
galaxies within 1 of the X-ray position with r — 3.6|J,m > 3 
and (independently) r > 22 and 3.6|J,m — 4.5|J,m> —0.1 was 
computed. For 3.6^im sources that are undetected in r (ap- 
proximately r > 25.9) the computed colour is a lower limit 
on the true colour. These values were compared to those 
computed for the z < 0.8 clusters and 1000 randomly placed 
apertures over the survey area. Figure |6] displays these val- 
ues for all CI and C2 clusters and cluster candidates in the 
9 deg^ area. As expected, the spectroscopically confirmed 
distant clusters are in the top-right of the distribution, along 
with many of the candidates. This is evidence that at least 
a subset of the candidates are at redshifts similar to, or 
even higher than, those of the spectroscopically confirmed 
clusters, with the candidate fields having a higher surface 
density of sources with the reddest 3.6^im — 4.5 ^im colours 
than the spectroscopically confirmed z > 0.8 clusters. 

The separation of cluster classes indicated in Figure [6] is 
supported by a two-sided Kolmogorov-Smirnov (KS) analy- 
sis of the surface density values (in this case using r — 3.6|J,m 
colours). A KS test between the z < 0.8 cluster sample and 
the 1000 random apertures result in a probability that the 
two samples are drawn from the same population of 0.74, 
which is sensible given the photometric criteria are designed 
to be sensitive to z > 0.8 galaxies. The corresponding KS 
probabilities between the z > 0.8 sample (confirmed and 
candidates) and the z < 0.8 clusters, the random apertures 
and the unknown sources are 5 x 10~^, 5 x 10~^ and 4 x 10~^ 
respectively - thus confirming the trend observed in Figure 

ii 

However, some overlap remains between the visually- 
assigned classes. For the low- versus high-redshift sources 
this is partly due to the fact that photometric uncertain- 
ties will result in low redshift galaxies exceeding the applied 
colour cut and vice versa - this blurring is expected to be 
most evident for clusters close to the effective redshift cut-off 
implied by the photometric criteria. In addition, while it is 
unlikely that any of the 57 z < 0.8 clusters are at ^ > 0.8 we 
do note that two of the z > 0.8 candidate clusters ultimately 
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Figure 6. The surface density of galaxies selected according to 
the photometric criteria described in the text and located within 
1' of the X-ray centroid. Sources are classified as 2; < 0.8 (blue 
squares), z > 0.8 spectroscopically confirmed clusters (solid red 
triangles), z > 0.8 candidate clusters (open red triangles) and 
unknown sources (green circles). 



result in photometric redshift estimates slightly less than 
z = 0.8. In general though, when considering the z > 0.8 
candidate clusters and the unknown extended sources the 
surface density analysis indicates that the optical-MIR data 
cannot provide an unambiguous assessment of these systems. 
Put another way, there is no straightforward threshold that 
can be applied to the surface density of optical-MIR selected 
galaxies along the line of sight to the extended X-ray source 
sample that will generate a complete sample of candidate 
distant cluster candidates with low contamination. We ad- 
dress this point further in Section |4] 



3.3 Identifying cluster red sequences 

In addition to computing the surface density of faint, red 
galaxies, colour magnitude diagrams (CMDs) and colour his- 
tograms in both r — 3.6^Lm and 3.6 ^im — 4.5^im were created 
for the confirmed and candidate z > 0.8 clusters and com- 
pared to the background colour distribution (e.g. see Figure 
|7|. Assuming that a cluster contains a significant number of 
passively-evolving galaxies, the presence of a red sequence 
should confirm a cluster identification and the characteris- 
tics of that sequence should indicate an estimated redshift 
for the cluster when compared to stellar population models 
(Figure |8|. The analysis indicates that the confirmed clus- 
ters, located at 0.8 < z < 1.2, display red sequence colours 
broadly consistent with the expectation of a simple model 
of an old, passively evolving stellar population considered at 
the confirmed spectroscopic redshift (see Figurejsjfor further 
details). A number of the candidate clusters at z > 0.8 also 
display red sequences consistent with the expectation of a 
high redshift passive stellar population. However, it is also 
clear that a subset of the candidate clusters do not display 
red sequence colours consistent with this simple model of 
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Figure 3. Colour images of spectroscopically confirmed clusters (see Table[T]). Images are composed of r2;3.6|Lim data. The white squares 
indicate spectroscopically confirmed galaxies at the cluster redshift. The white contours indicate X-ray emission. Images are 2' on a side 
with standard astronomical orientation. 



a high redshift stellar population - specifically they appear 
to be systematically bluer in r — 3.6^Lm than the expected 
colour for 1 < z < 2. This effect can be understood by 
noting (as indicated in Figure [7| that a number of sources 
detected at 3.6^Lm in the candidate clusters are undetected 
in CFHTLS Wl r-band data. The indicated r - 3.6M,m red 
sequence colours are therefore lower limits and will under- 
estimate the true colour. A further limitation is that the 
3.6|J.m— 4.5|J,m red sequence colour is largely degenerate with 
redshift at z > 1.3. The above steps indicate (often strongly) 
that each confirmed and candidate cluster is associated with 
an overdensity of galaxies consistent with z > 0.8. However, 
although the optical and MIR imaging data are effective at 
determining the presence of over-densities of high-redshift 
galaxies they only provide limited information on the red- 
shift associated with the red sequence location. 



3.4 Photometric redshift analysis 

Having searched for evidence of red sequences in colour mag- 
nitude diagrams (CMDs), we extend our analysis by remov- 
ing the assumption that any distant cluster galaxy popula- 
tion has a significant red sequence, and allow the possibility 
that the galaxies associated with a cluster can have a range 
of spectral energy distributions corresponding to a range of 
star formation histories. 



We use the public code Le Phare ( [Arnouts et al.|[2002 
Ilbert et al.|2006 ) to estimate the photometric redshifts. Le 
Phare is based on a standard template fitting procedure. 
The templates are redshifted and integrated through the 
appropriate transmission curves. The photometric redshifts 
are obtained by comparing the modelled fluxes and the ob- 
served fluxes with a merit function. We run the code us- 
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Figure 4. Colour images of candidate clusters (see Table [T]). Images are composed of rzS.6\iin data. The white contours indicate X-ray 
emission. Images are approximately 2^ on a side with standard astronomical orientation. 
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Figure 4 - continued Colour images of candidate clusters (see Table [TJ. Images are composed of ^ 
indicate X-ray emission. Images are 2' on a side with standard astronomical orientation. 



2[3.6|Lim] data. The white contours 
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Figure 5. Colour images of marginal or unknown extended X-ray sources (see Table [TJ. Images are composed of r2;3.6M.m data. The 
white contours indicate X-ray emission. Images are 2' on a side with standard astronomical orientation. 
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Figure 7. Optical-MIR photometry for cluster candidate 21. All 
sources within V of the X-ray position are indicated. Top panel: 
red points indicate sources undetected in r. Bottom panel: the 
histograms indicate sources detected in r (blue), undetected in r 
(red) and total source counts (black). Each histogram is scaled to 
be visible in overlay. The dashed line shows the background colour 
distribution. This is formed by computing the colour histogram 
of all sources more than V away from each C1+C2 source in the 
sample. The histogram amplitude is then scaled by the ratio of 
the areas of the cluster and background areas. 
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Figure 8. The distribution of cluster red sequence locations in 
the r — 3.6|J.m and 3.6|J.m — 4.5|J.m plane. Blue points indicate 
spectroscopically confirmed clusters at 2; > 0.8 and red points 
indicate candidate clusters at z > 0.8. Errors are fixed at 0.1 
mag. in 3.6fxm — 4.5|xm and 0.5 mag. in r — 3.6|Lim. The black line 
indicates the locus followed by a passively evolving 1 Gyr solar 
metallicity burst of star formation occuring at 2;^ = 5 ( |BruzuaT| 
|fc Chariot]|2003[| . The colour of this model stellar population at 
redshifts 1, 1.5 and 2 are indicated respectively by the black circle, 
triangle and square. 



ing exactly the same configuration as used in the COSMOS 
fiel d (Ilbert et al. 2009). The set of templates was generated 
by Polletta et al. (2007) for the elliptical and spiral galax- 
ies. Twelve blue templates generated with |Bruzual Char-| 
lot (2003) were added. Four different dust extinction laws 
were applied fPrevot et aL||1984| [Calzetti et al.| '2000), and 
an additional bump at 2175A, depending on the considered 
template. Emission lines were added to the templates using 
relations between the UV continuum, the star formation rate 
and the emission line fluxes ( [Kennicutt ^ 1998 ) . Moreover, an 
automatic calibration of the zero-points was performed using 
a sample of 650 spectroscopic redshifts within the photomet- 
ric data area drawn from the XMM-LSS sample described 
by Adami et al.| ( j2011| ) . The calibration is obtained by com- 
paring the observed and modelled fluxes ( Ilbert et al.|2006' ) 
at known spectroscopic redshifts. 

Having demonstrated that the optical and MIR colours 
of candidate 2; > 0.8 galaxies provide only a relatively inac- 
curate estimate of the redshift of a given confirmed or can- 
didate cluster, the photometric redshift analysis described 
here is limited to those confirmed and candidate clusters 
with additional NIR photometry in at least the J and Kg 
bands (see Table [l] for a list of clusters with JK photome- 
try). At 2; > 0.8 the JK bands sample the rest frame optical 
spectral energy distribution (SED) including the prominent 
D4000 feature which provides a strong JK colour signature 
- and thus strong constraining power in a photometric red- 



shift analysis - as a function of redshift for galaxies com- 
posed of evolved stellar populations. Figure |9] displays the 
photometric redshift histograms for all clusters with avail- 
able JK data. For clusters with spectroscopic redshifts, data 
are plotted for all galaxies within 30^^ of the X-ray source 
that are brighter than the i^s-band completeness limit of 



each data set (see Section 2.2). For clusters with photomet- 
ric redshifts (see below) Zphot < 1-2 and Zphot > 1-2, data 
are plotted for all galaxies within 1' and 30^^ respectively of 
each X-ray source that are brighter than the corresponding 
Ks-hand completeness limit. 

We represent the redshift density function for each clus- 
ter using both a standard histogram (6z = 0.1) and a 
variable kernel density estimation (VKDE) approach. The 
VKDE approach represents the contribution of each galaxy 
to the redshift density function as a Gaussian of width 
a — 0.03 X (1 -\- z) and unit area. A first estimate of the 
photometric redshift of each cluster is determined by identi- 
fying visually the peak associated with the X-ray source in 
Figure [9] We then compute the cluster photometric redshift 
as the mean of the VKDE-weighted redshift distribution, i.e. 



Zcluster 



J z K{z) dz 
J K{z) dz 



(1) 



over the local interval where the VKDE distribution K{z) 
exceeds 0.5 x K{zpeak)- The resulting Zduster for all can- 
didate clusters is displayed in Table [l] For the five spec- 
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troscopically confirmed clusters with NIR data coverage we 
can compute the effective photometric redshift error of this 
method as 

. 1/2 

N ^ V (1 + Zspec) 



-y 



(2) 



which yields Gz — 0.025. This error will naturally increase 
as one extends this approach to the candidate clusters which 
typically represent less clear redshift peaks composed of 
smaller numbers of galaxies. Figure |9] also displays photo- 
metrically selected cluster galaxies in field of each cluster. 
Cluster galaxies are selected as occupying the local inter- 
val K{z) > 0.5 X K{zciuster)- Flgurc [o] furthcr displays the 
J — Ks CMD for photometrically selected members of each 
cluster. 



4 DISCUSSION 



Cluster red sequences: implications for galaxy 
assembly history 



4.1 



In Section [3^ we investigated the limitations of the r— 3.6^tm 
and 3.6^im — 4.5|am data in the determination of the loca- 
tion of cluster red sequences. We return to this issue with 
the NIR data and the results of the photometric redshift 
analysis in hand and investigate two approaches to compute 
the colour sequence of cluster member galaxies selected by 
a) photometric redshift and 2) statistical background sub- 
traction on the NIR colour magnitude plane. 

4.1.1 Computing the red sequence colour employing 
photometric redshifts 

Figure [9] displays the CMDs of cluster galaxies in confirmed 
and candidate clusters selected by photometric redshift (or 
spectroscopic redshift if available) and proximity to the X- 
ray source. In order to identify the location of any red se- 
quence in each cluster we apply the method outlined by |Fass-| 
bender et al. ( 2011| ) whereby one identifies the third reddest 



galaxy in each CMD (3RG) and selects galaxies displaying 
J — Ks colours within 3RG±0.3. We then compute the clus- 
ter "red sequence" colour as the median of this set of galaxies 
iCmed) and the spread as either the colour of the galaxies 
enclosing 68% of the distribution (i.e. Gc = 0.5 x [C84 — Cie]; 
for A/" > 3) or the full colour range (for N ^ ?>). We apply 
this method to all confirmed and candidate clusters with 
J — Ks imaging and display the resulting Cmed and Gc val- 
ues in Figure |10| For cluster 19 we set Cmed equal to the 
colour of the candidate BCG and Gc = 0.1. In many cases 
this approach identifies a viable, yet often poorly populated, 
red sequence. However, we note that this approach does not 
identify a clear red sequence for candidate clusters 16, 17 
and 18. We address this result in more detail in the Section 

For those clusters displaying a clear red sequence con- 
sistent with Cmed and Gc we compare the red sequence 
colour versus redshift data to a set of representative colour- 
redshift loci generated using the BC03 spectral synthesis 
code (Bruzual Chariot 2003, ). The plot indicates that 
the location of the putative red sequences of the confirmed 
and candidate clusters follow a clear locus in colour-redshift 



1.6 



I 



-T 1 1 1 1 1 1 I I 1 1- 




2.5 



redshift 



Figure 10. Cluster red sequence colour computed for spectro- 
scopically confirmed (solid squares) and candidate (open squares) 
clusters. The colour-redshift data are compared to model predic- 
tions generated using the BC03 spectral sysnthesis code (Bruzual] 
|fc Chariot |2003| ). Each model features either a 1 Gyr burst (solid 
line) or r = 1 Gyr (dashed line) solar metallicity stellar popu- 
lation with a standard Salpeter IMF. Models are displayed for 
formation redshifts of 10 (red), 5 (green), 4 (cyan) and 3 (blue) 
respectively. 



space that can be described by the simple passive evolution 
of an old stellar population. We do not attempt to fit the 
stellar population parameters best describing the data be- 
yond noting that solar metallicity models arising from a 1 
Gyr burst of star formation at > 5 appear to be favoured 
over either extended {r — 1 Gyr), sub-solar (Z = OAZq) or 
younger (z/ < 5) bursts of star formation. It is worth noting 
at this point that the strength of this conclusion rests upon 
the red sequence colours of the candidate clusters at z > 1.5. 
In one sense this demonstrates the important leverage that 
distant clusters place upon our knowledge of galaxy evolu- 
tion in dense environments. In another sense however, this 
result can only be viewed as tentative pending spectroscopic 
confirmation of these clusters. 



4.1.2 Computing the red sequence location employing 
background subtraction on the CMD plane 

An alternative approach to identifying the red sequence of 
candidate cluster galaxies is to investigate their distribution 
directly on the colour magnitude plane. The main issue is 
to isolate the potentially weak cluster signature from the 
"background" of non-cluster galaxies along the line of sight. 
To achieve this we bin the colour magnitude distribution 
of galaxies located within 1^ of each cluster X-ray centroid. 
Galaxies are binned on the J — Ks versus Ks plane with bin 
dimensions of 0.2 and 0.5 mag. respectively. The CMD of 
non-cluster galaxies (termed the model background here) is 
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Figure 9. Photometric redshifts and NIR photometry for each cluster with available JK data (see text for additional information). 
Left panel: photometric redshift histogram for all sources located close to each X-ray cluster centroid. The red curve indicates the 
VKDE distribution and the dotted line is the appropriately scaled background. The vertical red tick mark indicates the computed 
photometric redshift for each cluster. For spectroscopically confirmed clusters the vertical black tick mark indicates the spectroscopic 
redshift. Middle panel: iJK image of each cluster. Each image is 2' on a side and follows standard astronomical orientation. The white 
contours indicate the X-ray emission in each field and the green circles indicate photometrically selected cluster members. Right panel: 
colour magnitude diagram for sources within 1' of each cluster X-ray centroid. For spectroscopically confirmed clusters the black squares 
indicate spectroscopically selected cluster members. For candidate clusters the black squares indicate photometrically selected members. 
Points indicate all other sources. The horizontal red and dotted black lines indicate the computed red sequence location and associated 
uncertainty. The angled dotted line indicates the photometric completeness in each field. 



formed by selecting all sources at > 1^ from each cluster or 
cluster candidate and scaling the resulting distribution by 
the relative cluster and background sky areas. 

Although the resulting background subtracted CMD for 
each cluster displays an identifiable overdensity of faint, red 
galaxies, the signal is accompanied by variations in the true 
background associated with each cluster (both Poisson and 
cosmic in nature). Upon subtraction of the scaled model 
background these variations persist as residual positive and 
negative signatures. However, their statistical distribution 
should average to zero over all cluster fields and, in an 
attempt to reduce their impact, we stack the binned sub- 



tracted CMDs to investigate the average cluster CMD. Prior 
to stacking we shift each cluster distribution on the colour 
magnitude plane to account for the /c-correction and dis- 
tance dimming of an evolving stellar population at each clus- 
ter redshift. We correct each cluster to a common redshift 
applying an apparent colour and magnitude shift based upon 
the evolution of a 1 Gyr solar metallicity burst of star forma- 
tion occuring at = 10 and described by a Salpeter IMF. 
Using this same model we have confirmed that the colour 
terms between the different NIR filter systems employed are 
small compared to the photometric zero point errors for each 
cluster. We stack all cluster and cluster canditates as follows: 
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Figure 9 - continued 
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spectroscopically confirmed clusters at z > 0.8 (5 systems), 
cluster candidates at 0.8 < Zphot < 1-2 (4), cluster candi- 
dates at Zphot > 1-2 either with a clear red sequence (4) 
or without (3, i.e. clusters 16, 17 and 18). We also stack 
the CMD of six control fields constructed in order to test 
the null hypothesis that each cluster candidate is false. We 
take the location of six clusters observed with HAWK-I and 
shift the cluster centroid to an adjacent detector (approx- 
imately a 3^ shift). We then repeat the stacking procedure 
using these new centroids and apply the same colour and 
magnitude shifts applied to galaxies selected according the 
original cluster locations. 

Figure [TT] displays the colour histrogram generated from 
each stack by summing the CMD along the magnitude axis 
with the restriction Kg < 22.5, for J — Kg < 0.5, or 
Ks < 23.0 — {J — Ks) otherwise, to consider the photo- 
metrically complete region of the CMD. Figure pT] confirms 
that there is little difference in the intrinsic red sequence 
properties between spectroscopically confirmed clusters at 
z ^ 1 and the candidate clusters. More importantly, the sub- 
set of candidate clusters lacking apparent red sequences on 
the basis of photometric redshift selection have been shown 
to have average red sequences statistically identical to the 
remainder of the candidate sample following the stacking 
analysis. We speculate that these clusters display a range 
of star formation histories that are not well described by 
the available SED templates used in the photometric red- 
shift analysis. The resulting photometric redshift peak will 
be broadened by this systematic uncertainty and attempt- 
ing to select cluster galaxies using the photometric redshift 
method described in this paper will result in a greater level 
of background contamination relative to clusters with well 
modeled SEDs. Thus the already weak red sequence may be 
diluted further by the increased effective background along 
these sight lines. This effect can only be verified once spec- 
troscopic redshifts are available for these clusters. All of the 
stacked cluster colour distributions are clearly real and sig- 
nificant when compared to the null distribution. For the 
clusters without an individual red sequence the fact that 



only 3 systems contribute to the average supports the as- 
sertion that the majority (if not all) are real as otherwise 
the signal observed in the average histogram would be very 
much diluted. We therefore conclude that each of the clusters 
presented in this sample is real in that the extended X-ray 
source is associated with a galaxy population clustered both 
spatially and in colour, whose spectra, photometric redshifts 
and/or colours are consistent with z > 0.8. 

4.2 The cluster red fraction 

The cluster red fraction provides a simple statement of the 
population mix of cluster galaxy members selected by colour 
and magnitude. It is conceptually identical to the cluster 
blue fraction computed by Butcher and Oemler (1984) yet 
here we focus on the the red galaxy component in order to 
highlight the contribution of the cluster red sequence to each 
cluster in our sample. We compute the red fraction as the 
ratio 

Nr ^cluster 

-Nr, 

back /n\ 

JR ITr ITr (^) 

^^T, cluster — ^^T,back 

where Nr denotes the number of galaxies satisfying J—Ks> 
3RG — 0.3 and Ks brighter than a evolving model early- 
type galaxy. The model assumes a galaxy of Kg = 22.5 at 
z = 1.5 described by a IGyr solar metallicity burst of star 
formation occurring at Zf — 10 and described by a Salpeter 
IMF. The corresponding Kg limit at the redshift of each 
cluster is computed accordingly. Nt denotes the total num- 
ber of galaxies satisfying the appropriate magnitude limit. 
Nciuster includes all galaxies within V of the X-ray centroid 
and Nhack includes all galaxies at greater than V from the 
cluster centroid with the value scaled to match the relative 
areas of the cluster and background samples. 

The red fraction values for all confirmed and candidate 
clusters with J — Ks photometry are displayed in Figure 
[12] The figure indicates that the cluster sample displays a 
range of red fraction values ranging from clusters almost 
wholly dominated by the red sequence (/h ^ 1) to those 
with low values, i.e. Jr < 0.2. The observation of a wide 



© 2011 RAS, MNRAS 000,[Tp5l 



Distant galaxy dusters in the XMM-LSS survey 21 



-I — I — I — I — I — 1- 







1.2 I 1 — r-n — \ 1 1 1 — — p- 



o.e - 



I 13 1 



0.6 - 



0.4 



0.2 



-1 \ \ 1- 



"I E] 1 



I = 3 1 



_j I I 



0.5 



1 

J-K, 



L5 



1.5 



redshift 



2.5 



Figure 11. Colour histograms generated by summing over the 
stacked CMD histograms for each cluster subsample: spectroscop- 
ically confirmed clusters (solid black; 5 systems), candidate clus- 
ters with Zphot < 1-2 (green; 4 systems), candidate clusters at 
Zphot > 1-2 with evidence for individual red sequences (red; 4 sys- 
tems) and without (blue; 3 systems), 6 random locations within 
the HAWK-I fields (dashed black). The numbers represent the av- 
erage number of galaxies at each colour brighter than the applied 
brightness cut. See text for further details. 



range of red fraction values supports the assertion that the 
compilation of a complete sample of X-ray selected distant 
clusters can provide a relatively unbiased view of galaxy 
populations in such systems. A comparable analysis of the 
red fractions in a non- X-ray selected distant cluster sample 
has not yet been performed - which is unfortunate as the 
results of such a study would provide a valuable perspective 
on the wavelength dependent biases affecting distant cluster 
identification. 

The range of red fraction values displayed by the dis- 
tant cluster sample is similar in extent to that of compa- 
rable mass clusters at redshifts z ~ 0.3 (Urq uhart et ah] 



2010| ) also studied within the XMM-LSS survey (note that 
we discuss mass estimation based upon X-ray flux measure- 



ments in Section 4.3). However, the XMM-LSS distant clus- 
ters will increase in mass by a factor ^ 5 between a redshift 
z = 1.5 and z < 0.3 ( |Boylan-Kolchin et al.||2009| . Such 



low redshift clusters of mass > 5 x 10 M© typically display 
dominant, bright red sequence populations with /r > 0.8. 
If one assumes that the evolution of galaxies onto the red 
sequence occurs following the rapid cessation of star forma- 
tion (quenching) - for example as a result of ram pressure 



stripping (Gunn and Gott 1972) or galaxy-galaxy interac- 
tions ( Dressier et al.||1994 ) - then the large observed range 
of red fraction values displayed by the XMM-LSS distant 
cluster is consistent with the scenario whereby the galaxy 
populations have been caught in a variety of states trans- 
forming between active star forming environments (where 



Figure 12. The red fraction for all confirmed (solid squares) and 
candidate (open squares) z > 0.8 clusters with J— Kg photometry. 
Error bars on the red fraction are Poissonian. 



the red fraction is low and comparable to that of the field) 
to a "red and dead" environment typified by 2; < 1 massive 
clusters. 



4.3 X-ray fluxes and cluster masses 

A key motivation for identifying the most distant clusters 
is to determine their global properties and thereby reveal 
the details of cluster evolution. One of the most important 
properties is the total cluster mass. Unfortunately, this is 
not observable directly but it can, with certain assumptions, 
be determined from other observables including the X-ray 
flux (luminosity) of a system. While one could employ well- 
studied scaling relations at 2; < 1 in order to convert a flux 
measurement to one of mass, the extrapolation of these scal- 
ing relations to 2; > 1 is fraught with uncertainty. We there- 
fore adopt an alternative approach whereby we compare the 
X-ray flux observed from clusters of known redshift (either 
spectroscopic or photometric) to the flux expected from a 
model cluster of specified properites. This approach offers 
the reassurance that the model assumptions are defined and 
compared to the two observables (flux and redshift) in as 
clear a manner as possible. 

Figure ^] compares the flux values for individual clus- 
ters to the flux of model clusters computed as a func- 
tion of cluster mass and redshift. The model assumes the 



luminosity-temperature relation of Arnaud & Evrard ( 1999 ) 



with self-similar evolution. The mass-temperature relation 



is taken from Arnaud, Pointecouteau & Pratt (2005) with 
6 = 200 and massive clusters (extrapolated down to the en- 
tire mass range) with self-similar evolution. A comparison of 
the flux values for the individual clusters to the model indi- 
cates that the clusters display an approximate mass limit of 
6 X 10^^ to 1 X 1O^^M0. All clusters display masses inferred 
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Figure 13. X-ray flux measurements for clusters with conflrmed 
spectroscopic redshifts z > 0.8 (flUed triangles) and candidate 
clusters with photometric redshifts Zphot > 0-8 (open squares). 
The contours indicate the expected flux versus redshift properties 
of model clusters of given mass (see text for details). 



from this comparison less than 4 x lO^^M©. We refer to this 
as the baseline model in the following text. 

It is clear that adopting a different set of model assump- 
tions will affect the mass estimates returned by the above 
analysis. Figure[l4]examines the extent to which adopting an 
alternative set of scaling relations influences the estimated 
cluster mass. We retain the assumption that scaling laws 
evolve in a self-similar manner and compare two alternative 
approaches to our baseline model described above (Figure 



14 



black solid lines): 1) a model which replaces the M200 — T 
relation for a M500 — T relation taken from Sun et al. (2009), 
valid down to 1 keV, and assuming the L500 — T relation 
described by Pratt et al. (2009), see Figure 14 (blue dot- 



ted lines). Model 2) considers the flux-redshift relation for 
self-similar clusters following the scaling laws described in 
Vikhlinin et al. 



(2009) - see Figure 14 (red dashed lines). 
Each of these alternative models generates mass estimates 
for clusters of given flux and redshift which are a generally 
a factor 2 lower than those generated applying the baseline 
model. Figure [15] examines the extent to which allowing a 
given scaling relation to evolve with redshift affects the esti- 
mated cluster mass. Once again, the black lines in Figure 
indicate the baseline model described by self-similar evolu- 
tion. Alternative evolution prescriptions include that of Re- 
ichert et al. (2011; blue dots), Clerc et al. (2012; red dashed) 
and the evolving scaling relations derived by Vikhlinin et al. 
(2009; green dot-dash). In this case, introducing alternative 
assumptions regarding the evolution of scaling relations can 
generate mass estimates for clusters of given flux and red- 
shift which are up to a factor 2-4 greater than those gener- 
ated applying the baseline model. The above considerations 
indicate that models which consider potential variations in 
cluster scaling relations and their evolution can generate 
mass estimates for z > 0.8 clusters which vary by up to an 
order of magnitude. However, given the shape of the con- 




Figure 14. A comparison of the predicted flux versus redshift 
behaviour of model clusters of varying assumed scaling relations. 
Black solid lines indicate the baseline model of iArnaud EvrardI 
(1999) and| Arnaud, Pointecouteau Pratt] ( |2005t . Blue dotted 
lines indicate a model assuming a M500 — T relation from Sun et 
al. (2009) and a L500 -T relation from Pratt et al. (2009). The red 
dashed lines indicate a model employing scaling laws described in 
Vikhlinin et al. (2009). See text for further details. 




oT^ o 




Figure 15. A comparison of the predicted flux versus redshift be- 
haviour of model clusters of varying assumed evolutionary proper- 
ties. The black solid lines indicate the baseline model described in 
the text which assumes self-similar evolution. Alternative evolu- 
tion prescriptions include that of Reichert et al. (2011; blue dots), 
Clerc et al. (2012; red dashed) and the evolving scaling relations 
derived by Vikhlinin et al. (2009; green dot-dash). 
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tours displayed in Figures [T3j [T4j and [15] it is apparent that 
the relative masses of the z > O.S clusters are likely to be 
robust against uncertainties in the assumed cluster scaling 
relation model. 



4.4 The abundance of distant X-ray selected 
clusters 

Having compiled a sample of z > 0.8 confirmed and can- 
didate clusters it is instructive to compare their observed 
abundance to the number predicted using a calculation 
based upon all relevant cosmological, cluster scaling relation 
and selection considerations. We computed the expected 
abundance of distant clusters by assuming a [Tinker et al.| 
(2008 ) redshift-dependent mass function. These authors pro- 
vide a functional form of the halo mass distribution cali- 
brated using numerical simulations up to z ~ 2.5. For this 
calculation, we assume ag — 0.787 (D unkley et"al]|2009| ). 



Halo masses are then converted from M2006 to M200C follow- 



ing Hu & Kravtsov ( 2003 ) assuming a NFW mass profile and 



a concentration model from Bullock et al. (2001). The con- 



version from masses to X-ray observables is performed based 
upon the same baseline M2ooc—Tx and Lx —Tx scaling laws 
as quoted above. We assume a scatter of (T\^Lx\Tx — ^-6 in 
the Lx — Tx relation. Both scaling laws are assumed to 
evolve self-similarly. The XMM-LSS C1+C2 selection func- 
tion (Pacaud et al. 2006) accounts for the probability of 
detecting a cluster given its X-ray observables, namely its 
[0.5-2] keV count-rate and its apparent core radius (corre- 
sponding to a /3 model with /3 = 2/3). For this purpose, 
model count-rates are estimated assuming an APEC spec- 
tral model with abundance Z — O.3Z0. The core radius is 
assumed to scale with the halo radius following a simple 
scaling relation ( Clerc et al.|2011 ): i^c = 0.24 x R^qq. 

The final expected redshift distribution is integrated 
in various redshift bins: we predict 28.5 clusters lying at 
z > 0.8, 5.2 at ^ > 1.3 and 2.5 at z > 1.5 in 9 deg^ of 
surveyed area. By comparison, the sample of high redshift 
clusters presented in this work consists of up to 22, 7 and 6 
clusters at z > 0.8, 1.3 and 1.5 respectively, thus showing a 
rough agreement with model predictions. However we note 
a slight 2a) excess of 2; > 1.5 clusters in our sample. 
If confirmed, this excess could be due either to an overes- 
timate of photometric redshifts or to a unaccounted selec- 
tion bias that would arise due to e.g. an increased contam- 
ination of the X-ray flux from unresolved AGNs in these 
objects. Furthermore, model uncertainties also impact the 
predicted number of high-redshift clusters. In particular, re- 
laxing the self-similarity constraint on the evolution of the 
cluster mass-luminosity relation can lead to considerably dif- 
ferent predictions ( [Pacaud et al.||2007| ). 

A further comparison of interest is that between the 
surface density of XMM-LSS distant clusters and that of 
other X-ray selected distant cluster samples present in the 
literature. The surface density of z > 0.8 clusters in the 
XMM-LSS survey is 2.3 =b 0.5deg~^ above a nominal mass 
limit of - 1 X lO^^M©. This may be compared to a figure 
of 15 z > 1.1 clust ers in approximately 1 deg^ reported by 
Bielby et al.|(|2010l above a mass limit of 1 x lO^^M© In 



(though specific details regarding the selection function are 
currently unavailable). The surface density of clusters re- 
ported in jFassbender et al. (2011) lies between 0.3 and 1.3 
deg~^ with the exact value depending upon the subset of 



XMM observations analysed. We note that the 'Fassbe nder| 
et al. (^11) compilation containes two previoulsy published 



XMM-LSS clusters with the result that this comparison is 
largely but not completely independent. 

The variance observed in these reported surface densi- 
ties arises from differences in the techniques applied to select 
extended X-ray sources, confirm galaxy overdensities and 
subsequently compute photometric redshifts or to compile 
spectroscopic redshifts. A further point worth noting is that 
distant clusters are often detected in survey data originally 
compiled to study galaxy clusters at z < 1. As such, they 
represent a subset of marginal detections and whether they 
are subsequently classified and confirmed as distant clusters 
is a very sensitive function of the set of selection tests ap- 
plied to the data. With this in mind we have attempted to 
generate a complete sample of distant X-ray clusters in a 
manner that depends solely upon the X-ray data by per- 
forming an analysis of all extended sources in a subset of 
the XMM-LSS area. One cannot completely escape the re- 
quirement of input from other wavebands, e.g. the optical, 
NIR and MIR data employed in this paper. However, we 
have attempted as far as possible to select galaxy overdensi- 
ties associated with each extended X-ray source in a manner 
which is insensitive to the assumed star formation history 
of individual galaxies. 

More detailed follow-up of individual clusters, includ- 
ing spectroscopic and deeper X-ray imaging observations are 
currently underway. Ultimately combining the precise red- 
shift measurements into a self-consistent analysis including 
selection effects and model uncertainties (^Pacaud et al. 2007 



addition Fassbender et al.| ( ^2011) present a compilation of 22 
0.9 < z < 1.6 clusters detected in up to 79 deg^ of archival 
XMM observations with a mass limit of 1 - 2 X lO^^M© 



Maughan et al.|2012|[Reichert et al.|201l] [Clerc et al.|201lf , 
we will be able to draw firm conclusions from the abundance 
of high-redshift clusters. 



5 CONCLUSIONS 

The analysis presented in this paper effectively completes 
the assessment of 88 extended CI and C2 sources from ap- 
proximately 9 deg^ of XMM-LSS data. Of these sources 
59 display spectroscopic or photometric redshifts z < 0.8, 
21 sources display spectroscopic or photometric redshifts 
z > 0.8 (or colours consistent with the same redshift limit in 
the case of candidate clusters 23 and 24), and a remaining 
8 sources appear to be consistent with misclassified point 
sources or marginal detections. The sample also contains 
cluster 19 at a photometric redshift z = 1.67 which is in- 
cluded in this paper having been flagged as a potential dis- 
tant cluster at an earlier stage. The distant cluster sample 
is complete in that it represents (with the low-redshift and 
marginal source identifications) a complete account of a 9 
deg^ area of the XMM-LSS survey. This sample is generated 
from the X-ray data employing a quantitative selection func- 
tion ( Pacaud et al.|2()06 ) and it therefore permits a number 
of important questions in cosmology and galaxy evolution 
to be investigated. 

The complete nature of this sample is dependent upon 
the primacy of the applied X-ray selection procedures. Al- 
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though it is difficult to conceive of a targeted X-ray survey 
for distant clusters that does not employ information at ad- 
ditional wavebands (e.g. optical, MIR, etc.), the application 
of a simple selection threshold to identify high-significance 
clusters in these additional wavebands (e.g. the surface den- 
sity of optical-MIR selected sources) will either lead to an 
incomplete X-ray sample or a large rate of contamination 
(from low-redshift or spurious sources), c.f. Figure |6] These 
comments do not undermine the nature of optical-MIR se- 
lected samples of distant clusters - which are complete in 
terms of the applied selection criteria. Instead they reflect 
the fact that distant X-ray detected clusters must be con- 
firmed at other wavebands. In doing so with this paper we 
have attempted to perform as comprehensive an assessment 
as possible of each confirmed and candidate system with the 
aim of compiling a complete sample of X-ray selected distant 
clusters. 

It is important to recognise that the analysis presented 
in this paper represents only one stage in the creation a com- 
plete sample of distant X-ray clusters. Clearly, much of the 
interpretation as to the nature of each cluster rests upon the 
photometric redshift analysis and spectroscopic confirma- 
tion of the redshifts of these clusters must be considered as 
the next, important step. Furthermore, the extent to which 
point source emission from unresolved AGN (both within 
each cluster and superposed along the line of sight) modi- 
fies the appearance of a sample of distant X-ray clusters is 
not well understood. We intend to employ the superior an- 
gular resolution of the Chandra observatory to characterise 
the point source contribution to a representative sub-sample 
of the distant clusters presented in this paper. Only when 
these steps are complete will we have a better understand- 
ing of what constitutes a complete, robust sample of distant 
clusters. This XMM-LSS distant cluster sample represents 
an important resource and will form the basis for studies in 
both the growth of large scale structure and the evolution 
of cluster galaxies to be presented in forthcoming papers. 
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